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ABSTRACT: Different search strategies for supersymmetry have been employed by the 
LHC general-purpose experiments using early data. As proven by their early results, these 
strategies are promising, but raise the question of how well they will generalize for the 
future. We address this question by studying two thousand phenomenological minimal 
supersymmetric standard model parameter space points that come from a fit to indirect 
and cosmological data. We examine the 5<r discoverability of the points employing a typical 
ATLAS-type search based on missing transverse momentum (MET), a search based on an 
optimised Mt2 cut an d a combination of the two, taking into account standard model 
backgrounds. The discovery reach of the strategies can depend strongly on the systematic 
uncertainty in the background, subject to the stringency of the cuts and the details of 
the background simulation. By combining the MET and M?2 based strategies, with an 
integrated luminosity of 1 fb" 1 (10 fb" 1 ) at 7 TeV, 4-8% (42%) of the points are discoverable, 
depending on the systematic uncertainty on the background. At 14 TeV and with 10 fb _1 , 
96% of the points are discoverable. While the majority of points can be discovered by both 
strategies at y/s = 14 TeV and with 1 fb _1 , there are some that are left undiscovered 
by a MET search strategy, but which are discovered by the Mt2 strategy, and vice versa, 
therefore it is essential that one performs both in parallel. We discuss some of the factors 
that can make points more difficult to observe. 



KEYWORDS: Supersymmetry Phenomenology, Hadron Collider 



ArXiv ePrint: 1105.1024 



Contents 

1 Introduction 1 

2 Points investigated 3 

3 MET based search strategy 4 

4 Mt2 based search strategy 5 

5 Technical details 7 

6 f £di-dependent discoverability 9 

7 Exclusion limits from ATLAS 11 

8 Including systematic uncertainties 13 

9 Optimal M^f 14 

10 Difficult points 14 

11 Discussion 18 

12 Conclusions 20 



1 Introduction 

The search strategies for supersymmetry at the LHC experiments ATLAS [1] and CMS 
[2] have been designed to be as general as possible. At the time of their design they were 
most extensively tested by the experiments in simulations on 7 and 14 benchmark points, 
respectively. These points are thought to be representative for whole extended regions in 
the minimal supersymmetric standard model (MSSM) parameter space, but in fact belong 
to one particular constrained model called mSUGRA/CMSSM. Therefore it is possible 
that these points do not cover all of the MSSM parameter space satisfactorily. The search 
strategies (for example cuts) have been optimized to increase the efficiency of the searches 
in these particular cases. It is important to understand to what extent these searches will 
be able to discover SUSY in other, less constrained regions of the MSSM parameter space. 

The searches are based on the observation that in a wide range of SUSY models, squarks 
and gluinos are expected to cascade decay to jets, leptons or photons. The final state is 
expected to contain stable, weakly interacting (and therefore invisible) heavy particles, 
which make a good candidate for dark matter. These will escape detection, giving rise 
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to missing transverse momentum in the detector. SUSY searches at the LHC [3], [4] are 
heavily based on this signature, requiring large missing transverse momentum in order to 
discriminate standard model processes from supersymmetric events. In our paper, we call 
such a search a MET-based strategy. 

Another strategy examined in this paper is incited by the fact that many well-motivated 
models of new physics have a Z2 parity (for example i?-parity in super symmetry), which 
predicts that supersymmetric particles must be produced in pairs. Each of the two super- 
symmetric particles cascade decays within the detector to the lightest parity-odd particle, 
which is stable and assumed to be weakly interacting. i?-parity therefore gives rise to very 
particular event topologies. The Mt2 variable [5, 6] is a kinematic variable that uses ex- 
perimentally measurable quantities (the missing transverse momentum being one of them) 
efficiently, by exploiting the topology of the event, to extract information about the masses 
of the produced particles. Search strategies using Mt2 as a discriminating variable have 
been designed [7] and recently applied by the ATLAS experiment in [3]. 

The latest results from the ATLAS [3] and CMS [4] experiments have proven that their 
search strategies for supersymmetry are promising. However, we would like to understand 
how well these generalize to the future. Two recent studies [8, 9] have examined ATLAS 
MET based search strategies in a more general context by simulating signals for ~ 71k 
model points within the phenomenological MSSM (pMSSM). The landscape of the pMSSM 
parameter space has also been explored in a study by AbdusSalam et al. [10], who have 
performed a global Bayesian fit of the pMSSM to current indirect collider and dark matter 
data. The parameters of the pMSSM were constrained in the fit, resulting in inferences on 
sparticle masses. In the present paper, we take a sample of points from the fit of reference 
[10] with a probability proportional to their posterior probability function. Extending the 
work of [8, 9], which do not examine the discovery reach of M?2 based searches, in this 
paper we will use the points sampled from reference [10] to: 

• determine how likely not only MET-based searches, but also Mj^-based searches are 
to discover SUSY. 

• understand the successes and failures of the searches. 

• optimise the searches for the more general pMSSM scenario. 

The plan of this report is as follows. First we give a short account of the properties 
of the SUSY points examined, sampled by [10]. Then, the MET and Mt2 based search 
strategies and the corresponding event selection are explained in sections 3 and 4. We 
briefly present in section 5 the technicalities of simulating the events for different model 
points, a basic detector simulation, a discussion of the Standard Model backgrounds and 
provide a definition of discoverability. Section 6 shows how the SUSY discoverability varies 
with the choice of the M?2 cut for different centre of mass energies (7 and 14 TeV) and 
a variety of integrated luminosities, in the context of a search strategy based on Mt2 
alone and a combined search strategy based on both MET and Mt2- We discuss the most 
recent ATLAS exclusion limits on sparticle masses and compare them with the expectations 
from our study in section 7. The effect of systematic uncertainties in the background on 
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discoverability is detailed in section 8. Optimal cuts are proposed in section 9. We discuss 
what makes a point difficult to discover in section 10. Finally, section 11 contrasts our 
work on the Mt2 search strategy with the analysis from references [8, 9] performed on 
other ATLAS MET-based search strategies. We conclude in section 12. 

2 Points investigated 

A total of 1940 points were provided by the authors of [10]. They are sampled from the 
posterior probability distribution obtained through a Bayesian fit of electroweak, search and 
flavour physics data to the pMSSM. The pMSSM contains the most relevant 25 weak-scale 
MSSM parameters. The priors used in the analysis are flat in mass parameters except for 
the scalar mass parameters, which are flat in the logarithm. Inferences are made about 
these points without assuming a restrictive high-scale supersymmetry breaking model. The 



Table 1: Allowed ranges for the pMSSM parameters, tan/3 is the ratio of the MSSM Higgs 
vacuum expectation values. 

ranges of allowed values for some of the parameters are summarized in table 1. Five of 
the 25 pMSSM parameters varied were associated with the Standard Model: the top mass, 
the mass of the Z boson, the bottom quark mass and the strong coupling constant. The 
relic density of dark matter inferred from cosmological measurements proved to be the most 
important constraint. Various observables involving bottom quarks were used to constrain 
the model: the branching ratios of B$ — > B u — > tv and B — > X s ^f as well as the 

B a — B s mass difference and the B — > K*^ isospin asymmetry. The electroweak observables 
used were: the W boson mass, the mass and decay width of the Z boson, the effective weak 
mixing angle, Z— pole asymmetry parameters and the muon anomalous magnetic moment. 
Direct pre-LHC searches for sparticles and Higgs bosons were also used. The 25 dimensional 
parameter space was successfully scanned with the MultiNest algorithm [11, 12], achieving 
a statistically convergent fit. The resulting fit was far from being robust, since there were 
not enough data to constrain so many parameters, and this manifests itself in the high 
degree of prior dependence in the results. In the present paper, we are not concerned with 
the issue of fit robustness: we merely want to use a subset of the points that fit current 
data well in order to investigate the properties of proposed search strategies in a general 
context, as explained in Section 1. 

Together with the SPSla point [13], these points form the pool of scenarios studied 
in this report. A flavour of their mass spectra is given in figure 1. In the sampling from 
the pMSSM, the physical masses of the sparticles was constrained to lay below 4 TeV. In 



Parameter 



Allowed range 



Scalar masses 
Trilinear scalar coupling 
Gaugino masses 
tan /3 



[100 GeV, 4 TeV] 
[-8 TeV, 8 TeV] 
[-4 TeV, 4 TeV] 
[2,60] 
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Mass of lightest squark (GeV) Mass of x, (GeV) 



Figure 1: A flavour of the mass spectra of the 1941 SUSY points investigated. The * 
denotes the SPSla point. A * is missing from the first histogram since the heaviest squark of 
the SPSla point (m = 586 GeV) is much lighter than the corresponding squark in any of the 
other points examined. 

contrast, all seven mSUGRA/CMSSM benchmark points used by the ATLAS collaboration 
in [14] have a gluino mass below 1 TeV. Notice that the SPSla point, which is used widely 
as a benchmark point, is also in the low mass end of the spectra and is now excluded by 
the ATLAS experiment [3]. 

3 MET based search strategy 

The final state of supersymmetric events is expected to contain stable, weakly interacting 
(and therefore invisible) heavy particles. These will escape detection, giving rise to miss- 
ing transverse momentum in the detector. In this paper we study one possible discovery 
channel, in which we require 4 or more jets and leptons in the final state. The lepton 
requirement suppresses the (semi)leptonic decays of W, Z and ti, whereas the large number 
of jets suppresses the QCD background. The concrete selection for this channel follows the 
recipe laid out in [1]. We require: 

• at least four jets with transverse momentum pj- > 50 GeV, at least one jet must 
have p T > 100 GeV 

• missing transverse momentum p > 100 GeV 
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no isolated leptons with px > 20 GeV 

the smallest azimuthal separation between the jet direction and the missing transverse 
momentum, A<j) l]et i ,p T ) > 0.2 for the first three jets 



• effective mass M cS > 800 GeV 

• f T > 0.2 M eff 

• transverse sphericity St > 0.2 

The missing transverse momentum p is calculated as the negative sum of the trans- 
verse momenta of all visible objects within detector acceptance: 

i 

where we define the acceptance as having pseudorapidity \rj\ < 5 and only take into account 
objects with px > 0.5 GeV. a 

The effective mass M e ff is defined as the scalar sum of transverse momenta of the most 
energetic objects in the event - in this case four jets, together with the missing transverse 
momentum: 

71=4 

M cff = ^| PT «|+^ T (3.2) 
1=1 

The transverse sphericity St is defined as: 

* s ira (3 ' 3 > 

with Ai and A2 being the eigenvalues of the 2x2 sphericity tensor Sij = Ylk (pt^ (pt^ 
The tensor is computed using all jets with \rj\ < 2.5 and pt > 20 GeV. 

A simplified version of this strategy, suitable for early data, was employed by ATLAS 
for three of the four signal regions defined in reference [3] . These follow the same philosophy 
as reference [1], but require a reduced jet multiplicity and slightly looser cuts to suit the 
limited statistics. The fourth signal region of reference [3] employed a different selection 
variable described in the following section. 
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4 Mt2 based search strategy 

The Mt2 variable is a generalization of the transverse mass to the case of decays of pairs 
of particles [5, 6]. Each of the two decays will contribute to the final state with one visible 
and one invisible object. The transverse momentum and Pt'' 2 '* of each of the visible 

objects can be measured in the detector. However, only the total transverse momentum of 
the invisible objects can be measured as missing transverse momentum p^. Nevertheless, 

a 7? = — In [tan , where 6 is the angle between the particle 3-momentum and the beam axis. 
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(i) 

for each of the decays labeled by i = 1, 2 one can compute the transverse mass rrij, under 
the assumption of a massless invisible object: 



f( Pr (^ T ( ! ))E2p r W ^W_ 2pT W.^ T W (4.1) 

where is a guess for the true, unknown missing transverse momentum ■ The 

variable Mj-2 is defined by: 

M T2 (pT (1) ,PT (2) ,y> T ) = min {max (m^m?)} ( 4 - 2 ) 

The minimization takes place over all values of the two undetectable particles' possible 
missing transverse momenta fijS 1 ' 2 ^ consistent with the constraint ^2$ T = $ T - 

The variable Mt2 gives the greatest lower bound on the parent mass for a given event 
that is consistent with the observed momenta and hypothesised masses [15]. When plotted 
over all events with the same parent particles, one expects an endpoint to appear. For 
example, for ti production, we expect Mt2 to have an endpoint at the mass of the top quark. 
While the endpoint will be smeared by effects such as detector resolution, simulations [7] 
show that the vast majority of standard model background events have Mt2 < nh- For 
SUSY events however, where the actors are far heavier particles, one can expect Mt2 to 
have larger values. This allows good discrimination between signal and background for the 
region Mt2 > m top- h Motivated by this observation, the use of the Mt2 variable as a 
complement to the standard ATLAS search strategies has been advocated in [7] and is now 
being actively pursued by the ATLAS collaboration, for example in [3, 16]. 

Typical Mt2 spectra for simulated backgrounds and for the SPSla SUSY scenario 
are displayed in figure 2, reproducing results of [7]. These were obtained as described in 
section 5. The figure shows that for a Mt2 cut > 600 GeV placed too high, not only the 
background, but also the signal is suppressed. For values of the Mt2 cut too low, the 
background swamps the signal. Our goal is to strike a balance between these two extremes. 
Based on the SPSla point one can find the range of values for a cut on Mt2 which will render 
the SPSla model discoverable, i.e. which 1) provides a sufficient number of signal events 
after selection and 2) provides a small enough p - value for the background-only hypothesis. 
The meaning of "sufficient number" and "small enough" will be defined in the next section. 
SPSla is only one possible scenario of how nature could look, so when determining the 
optimal value for the cut on Mt2, we perform the same analysis on many other scenarios. 

ATLAS has implemented in one of the four signal regions in reference [3] an Mt2 
search strategy by requiring no leptons, a jet multiplicity > two, with the leading jet 
p T > 120 GeV, second jet p T > 40 GeV, f T > 100 GeV, A</> (jet, p^) > 0.4 and 
Mt2 > 300 GeV. The Mt2 variable is calculated using the two highest-p^ jets. The 
p T requirement makes an effective trigger, while the A(j) requirement reduces the QCD 
background. 

In this paper we take a simplified version of this analysis, with the intention of examin- 
ing the effect of a cut on Mj-2 alone. We only require two or more jets with p? > 50 GeV, 
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b We note that for fairly degenerate mass spectra, where the splitting between sparticles is less than the 
top mass, good discrimination may not be possible. 
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M T2 (GeV) 

Figure 2: Distribution of Mt2 for events with two or more jets with pr > 50 GeV (and no 
other cuts) for an integrated luminosity of 1 fb _1 at s/s = 7 TeV. 

compute Mt2 using the two highest-pr jets and then cut at different values of Mt2- We 
call the value of Mt2 at which the cut is set Mffy . No restriction is set on the number of 
leptons. 

5 Technical details 

The following operations have been performed in this analysis: 

• The spectrum information for each of the 1940 pMSSM points sampled by the authors 
of [10] and for the SPSla point was passed to Herwig++ 2.4.2 [17, 18] via the SUSY 
Les Houches Accord [19]. For each point we simulated 10000 inclusive SUSY events. 
A sample of background events was produced, containing 10 7 events of each tt, W — > 
Iv + jets, Z —7- l + l~ + jets, Z — > vv + jets. For QCD we produced 10 7 events for 
each slice of px of hard scatter, where the slices lie between 17, 35, 70, 140, 280, 560, 
1120 GeV, with the last slice starting at 2240 GeV. This is to ensure we have sufficient 
simulated events at high pt- The contribution from diboson+jets is expected to be 
very small according to [1]. For each of the SUSY scenarios and the background, we 
further performed the following: 

• For each event, we clustered hadrons with fiducial pseudorapidity (\r]\ < 5) and mo- 
mentum (pt > 0.5 GeV) into jets using the fastjet [20] implementation of the anti-fcy 
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algorithm [21]. The E combination scheme with R = 0.4 and pT mm = 10 GeV was 
used. 



• We simulated the effect of the detector in a simplified way by smearing the energy 
and momenta of the jets using a Gaussian probability function for the majority of 
(1 — e) events. The width of the Gaussian depends on the energy of the jet such that: 



where Ej is the unsmeared jet energy. This resolution is typical of general-purpose 
LHC detectors [1, 2]. The distribution also has a low energy tail to account for the 
possibility of badly mismeasured jets, therefore the remaining fraction e = 1% of the 
jets is smeared with a probability density: 



where r = E/Ej is the ratio of the smeared jet energy to the true jet energy. The 
constant c\ is chosen such that the resulting smearing function of a jet of given energy 
agrees with full simulation results as in [1] and C2 is a normalization constant. 

The missing transverse momentum is calculated from the negative vector sum of the 
jet momenta (after smearing) and other isolated particles present. 

In addition, for the MET and Mt2 search strategy individually, we determined the 
number of events passing the cuts defined in sections 3 and 4. For the M?2 based strategy 
we allowed for a variable cut on Mt2- We then compared the number of signal events 
with the number of background events to determine whether the point in question would 
be discoverable at a given luminosity. For a SUSY point to be discoverable, we require on 
one hand at least 10 SUSY events to pass the cuts. On the other hand, the signal should 
be well above background, i.e. the p - value for the background-only hypothesis should 
be low enough. We will convert the p - value into an equivalent significance Z, defined 
such that a Z standard deviation upward fluctuation of a Gaussian random variable would 
have an upper tail area equal to p. That is, Z = <1>~ 1 (1 — p), with $ being the cumulative 
distribution of the Standard Gaussian. We then require for a discovery that Z > 5. We can 
include the systematic uncertainty in the background by using the methods from reference 



Having decided for each SUSY point individually for which range of M^f it would be 
discoverable and whether the MET search strategy would discover the point, we would like 
to know the answers to the following questions: 

• What is the overall optimal value for Mj^l 

• How is the optimal value of MjH£ affected as we include systematic uncertainties in 
the background? 




(5.1) 




(5.2) 



[22]. 
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• How is the optimal value of affected as the experiments accumulate more lumi- 
nosity? 

• How does the optimized Mt2 search compare to the MET search when taking into 
account the above? 

We will address these questions in the following sections. 
6 J £(f£-dependent discoverability 

At the time of writing, the LHC is scheduled to run at y/s = 7 TeV, where it is planned to 
accumulate a total of a few fb _1 . Figure 3a shows the number of points discoverable with 
different luminosities at 7 TeV as a function of M™ . For each integrated luminosity, two 
lines of the same style were drawn. The lower one shows the number of points discoverable 
with the Mt2 based search strategy alone. The upper of the two lines with the same 
style shows discoverability for a combination of the Mt2 an d MET based strategies, where 
a point is considered to have been discovered if it was discovered by either the Mt2 ot 
the MET-based strategy. The flat sections of the upper lines give the number of points 
discoverable with the MET based search strategy alone. If no points were discovered with 
the MET based strategy, the M?2 only and combined strategy lines overlap. Only the 
histograms corresponding to an integrated luminosity of up to 1 fb^ 1 have been or will be 
experimentally explored with a center-of-mass energy of 7 TeV. 

The same study can be applied for the LHC at its design center of mass energy of 
14 TeV. Figure 3b shows the discoverability in this scenario. Due to the increased center of 
mass energy, high discoverability is achieved with far lower integrated luminosity. 

Notice in figure 2 that our background sample is statistics limited at the high Mt2 
end, giving rise to a ragged, discontinuous spectrum. This translates in figures 3a-3d into 
possibly fake peaks of discoverability at large Mt2 values. Therefore, one should view 
critically predictions of discoverability for y/s = 7 (14) TeV for values of Mj^ above 
520 GeV (780 GeV). 

Once the discoverability profile and an optimal M™ value is found through simula- 
tions for a given integrated luminosity, it is not necessarily obvious if the same profile or 
optimal M™ value will apply for different luminosities. In the limit of large numbers of 
events and no systematic uncertainties in the background, the significance varies as cx 
This implies that as more luminosity is accumulated, discoverability will increase at low 
values, where high background is a problem for the significance. At the high MjH£ 
end, discoverability is not necessarily limited by the significance as much as by the fact that 
not enough signal events will have been recorded. Higher integrated luminosity increases 
discoverability for very high Mj^ values as well, since the number of signal events is pro- 
portional to luminosity. Considering these two effects together, we expect the "peak" of 
the discoverability distribution with respect to Mj^ to broaden and to shift to higher Mt2 
values. This explains the very flat maximum observed at high luminosity in figures 3a and 
3b. One can choose a low Mt2 value which optimizes the search at low luminosity - the flat 



- 9 - 




— 1 ab" 1 


cn 




- 1 ab 1 


o 10 3 
a. 




- 100 fb" 1 




F:." ... 


- 1001b' 


10 fb" 1 


U) 




10 fb'' 


1 fb 1 


overable 




1 fb 1 


100 pb 1 
— 35 pb 1 




100 pb' 
-35 pb ' 




disc 




i: J 






Number oi 


T" 


, i , , , i , , , 1 . 1 . i . . . i . . J i . . . i 





1000 1200 1400 
M^ 1 value (GeV) 



(a) ^/s — 7 TeV. Systematic uncertainties in 
the background have been neglected. 



1000 1200 1400 
M^ 1 value (GeV) 



(b) = 14 TeV. Systematic uncertainties in 
the background have been neglected. 




7 TeV 


n 


- 1 ab" 1 


'5 10 3 
a. 


— 100 fb" 1 


> 
tfl 


10 fb" 1 


cn 


1 fb 1 


3 


100 pb" 1 


2 10 2 


— 35 pb" 1 


Number of discovi 
o 



1000 1200 1400 
M^' value (GeV) 



(c) y^s = 7 TeV. The systematic uncertainty 
in the background was assumed to be 50%. 
Notice no points can be discovered with 
35 pb^ 1 




1000 1200 1400 
M^' value (GeV) 



(d) y/s = 14 TeV. The systematic uncertainty 
in the background was assumed to be 50%. 



Figure 3: Number of SUSY points (out of 1941) discoverable at \fs = 7 or 14 TeV versus 
Mj^* (bins of 20 GeV) for the Mti based strategy (lower line of same style) and the 
combined strategy based on Mti and MET (upper line of same style). For %fs = 7 (14) TeV, 
predictions in the region above 520 (780) GeV might not be accurate due to the limited 
statistics in the background. 



maximum will ensure that the same value will provide a nearly optimal strategy at higher 
luminosity as well. 

We have studied discoverability by incrementing the luminosity by a factor of 10 each 
time. However, this does not buy nearly the same factor in discoverability. This is un- 
derstandable since on one hand Z oc \[C at the optimal value of M™*- On the other 
hand, if ^fs < 2mLSP) where ttjlsp is the mass of the lightest supersymmetric particle, then 
sparticles are simply not produced, regardless of the luminosity. 

From the corresponding discoverability plots it becomes apparent that, in the absence 
of uncertainties in the background, the optimized Mt2 based strategy is better suited than 
the MET based strategy for integrated luminosities up to several inverse femtobarns, after 
which the MET based strategy performs better. 
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In an attempt to explain this, we found that the MET based search is too stringent 
and does not allow any background events at all and only very few signal events through. 
Therefore there are plenty of SUSY points which have fewer than 10 signal events passing 
the cuts of the MET based strategy. For these points, the discoverability is limited by the 
number of events, which grows linearly with the luminosity. On the other hand, the Mt2 
based strategy operating at the value of M^f giving maximal discoverability is significance 
limited and Z grows oc 

At low luminosity, not enough signal events pass cuts for the MET based strategy to 
be efficient. With increasing luminosity, the MET based strategy starts to dominate due to 
the linear growth of the number of events as opposed to the oc y/Z. growth of significance 
for the Mt2 based strategy. At high luminosity, of more than a few inverse femtobarns, 
this effect is so evident that the combination of the M^2 an d MET based strategies only 
discovers very few points more than the MET based strategy 

Figure 4 makes clear that while both the Mt2 and MET based strategies are good at 
discovering a large sub-space of the pMSSM, at high luminosities above 1 fb _1 the MET 
based strategy explores the region of heavy gluino masses, m g > 2 TeV, more effectively. 
Whilst the two strategies have a large overlap of points they discover, there still are regions 
in the (m§ — mughtest <?) plane which are only explored by one or the other strategy. There 
are a variety of reasons for this, which include the inability to set a single cut threshold on 
Mt2 which is optimal for all points and the relative sensitivity of each strategy to different 
jet multiplicities, jet pt and decay topologies. As an example of such a reason, we note that 
the Mt2 based strategy was designed for the topology of pair- producing SUSY particles, 
whose direct decays are assumed to produce the hardest two jets in the event. These are 
the jets that ought to enter the calculation of Mt2- However, we found it is possible for 
one or both of the two hardest jets in certain events to be initial or final state radiation. 
These events will be more likely to fail the M?2 than the MET cuts. Similarly, points 
favoring event topologies where the hardest two jets are daughters of SUSY particles will 
be discovered with the Mt2 based strategy more easily. Since the the two strategies are 
predominantly accepting events of different topologies, a combined strategy tends to be 
more powerful than the individual strategies. 

We conclude this section by noting that so far the LHC SUSY searches [3, 4] have not 
been strongly affected by multiple in-time pp collisions. In the future, due to the increase 
in the instantaneous luminosity, the effects of such pile-up will have to be investigated in 
dedicated studies, to include the modified response of the detector. 

7 Exclusion limits from ATLAS 

Turning the problem around, we ask: could one have expected the LHC to exclude SUSY at 
95% confidence level (had it been one of the points in this study) , with the integrated lumi- 
nosity it achieved by the end of 2010, i.e. ~ 35 pb _1 and ignoring systematic uncertainties? 
The most stringent constraints up to now come from the ATLAS collaboration, which ex- 
cludes in reference [3] at the 95% confidence level (CL) points with gluino or squark masses 
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Figure 4: Points in the (m g — Slightest g) - plane, discovered by either Mt2 alone, or MET 
alone, or by both MET and Mt2, assuming ^/s = 14 TeV, a luminosity of 1 fb _1 and a fixed 
systematic uncertainty in the background of 10%. The Mt2 search operates at its optimum, 
which for this luminosity and systematic uncertainties is at M^= 660 GeV. 

up to 500 GeV in simplified models containing only squarks of the first two generations, a 
gluino octet and a massless neutralino. 

For an MjO£ value between 200 and 240 GeV (which is our optimal value at the given lu- 
minosity) the combined strategy from this paper (after an integrated luminosity of 35 pb" 1 ) 
would have provided evidence for discovery with p - value < 5% for 54 points, including 
the SPSla point. These points should have been discovered or excluded at 95% CL by 
the LHC experiments in references [3, 4]. Indeed, these 54 points have squark and gluino 
masses that fall within the exclusion limits of the simplified model of the ATLAS analy- 
sis [3]. The simplified model assumed m^i = 0, which is the value for which discovery would 
be most straightforward, and contain no other light supersymmetric particles. Whilst it is 
true that our 'discoverable' model points all have (m^ghtest qi m g) values which lie within 
the exclusion region of the ATLAS analysis, the inverse is not true; there are many points 
in the sample examined in this paper with squark and gluino masses that lie within the 
ATLAS simplified-model exclusion region which would not be discovered with 35 pb _1 - for 
a variety of different reasons. Some characteristic sparticle masses of the five points which 
need least luminosity for exclusion are shown in table 2. The characteristics of some of the 
points that are found to be most difficult to discover with 10 fb _1 at y/s = 14 TeV are 



- 12 - 



discussed in section 10. The points studied therein do not fall within the ATLAS exclusion 
region, but the reasons for which they are difficult to discover also apply to the elusive 
points lying within the ATLAS exclusion region. 



Point 


m-g (GeV) 


Mass of lightest squark 
from the first two gen- 
erations (GeV) 


Mass of lightest squark 
from the third genera- 
tion (GeV) 


m Xl (GeV) 


n65 


520 


395 (ul,ur,c l ,c r ) 


848 (6x) 


229 


nl88 


560 


430 (d R ,u L ,s R ,c L ) 


437 (6 X ) 


198 


SPSla 


604 


548 (d R ,u R ,s R ,c R ) 


401 (k) 


97 


p29 


597 


406 (u L ,c L ) 


1955 (ti) 


183 


p780 


546 


398 (u L ,c L ) 


306 (6i) 


199 



Table 2: SUSY points that would have been excluded or discovered at 5cr with 
/ Cdt = 35 pb" 1 at y/s = 7 TeV with a combination of the Mt2 and MET based search 
strategies as described in this paper. In two of the columns we show in parentheses the 
lightest squark flavour belonging to the first two or the third generation. For the case of the 
first two squark generations, we group together several squarks with near-degenerate masses 
if their mass difference is less than 1 GeV. 



8 Including systematic uncertainties 

Unavoidably, there will be both theoretical and experimental uncertainties in the back- 
ground. We model the uncertainties in the background by assuming a conservative system- 
atic uncertainty of 50%, which is comparable to that reported for signal regions B and D 
in the latest ATLAS SUSY search [3]. 

For the study of the effect of systematic uncertainties we assume the background to 
be a Poisson process with mean b, where b itself is a Gaussian-distributed quantity with a 
standard deviation of max {50% X b, 1 event I . For each SUSY point we calculate in bins of 
Myj 1 the p - value for the background-only hypothesis, which again can be translated into a 
significance as described in reference [22]. With the altered prescription for calculating the 
significance, the discoverability is changed to the one displayed in figures 3c and 3d. The 
intermediate regime of a systematic uncertainty in the background of 10% is also considered 
in this study. The effect of including the systematic uncertainties in the background is to 
push up the optimal value of ■ 

Assuming a constant fraction of 50% systematic error in the background, regardless 
of the value of M^, is a rather naive assumption. However, the detailed determination 
of the uncertainties in the backgrounds will require both complete detector simulation 
and supporting measurements from the LHC data, and so can only be performed by the 
experimental collaborations themselves. Our assumptions on the uncertainties should not 
affect the qualitative discussion we present in this paper, however. 
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9 Optimal Ml 



We will consider the optimal value of the Mt2 cu t to be the one which maximizes the 
number of discoverable points according to the two conditions outlined in section 5. These 
conditions take into account both statistical and systematic uncertainties. For the optimum 
determination we have calculated the number of discoverable points in bins of the Mtpf in 
sections 6 and 8 and found the bin with the most entries. 

There are two caveats linked to our determination of the optimum, neither of which is 
directly relevant for the experimental measurements. Firstly, our optimal Mtpf is based on 
the limited number of events that we have simulated in the tail of the Mt2 distribution for 
both signal and background. Secondly, the overall systematic uncertainty will vary across 
the space of the selection variables, whereas we have assumed it to be constant. What will 
be done in practice by the experiments will be to determine most of the backgrounds and 
their uncertainties from other measurements of LHC data. 

Table 3 summarizes the results for the optimal Mj^ value in the cases of no systematic 
uncertainty in the background, or an assumed 10% or 50% systematic uncertainty. The 
general trend is that the optimal value of MjJjj goes up with increasing luminosity, as 
expected from the discussion in section 6. In figure 3b we observe that the maximum of 
discoverability is very flat in M™, , such that choosing an M™* value which maximizes 
discoverability at say 100 pb _1 will be approximately optimal for higher luminosities as 
well. Considering the above, the following values look reasonable for the choice of the 
initial value for low luminosities, in the context of a combination of the Mti and 

MET based search strategies: for running at 7 TeV, 320 GeV; for running at 14 TeV, 
400 GeV. 

We also notice that with 50% uncertainty in the background, the M?2 based strategy 
barely adds any discoverability to the MET results. 

10 Difficult points 

Some difficult SUSY points evade detection unless they are tackled with a high energy LHC 
with high luminosity. Why is this the case? There are in principle two types of reasons: 

• The sparticle masses are very large with a consequently low SUSY production cross- 
section; 

• The mass spectrum and branching ratios lead to decay topologies / kinematics with 
low values of the variables which the search strategies cut on. 

Reasons of the first type need to be tackled with higher center-of-mass energies and 
higher luminosities. Reasons of the second type need to be tackled by improving the search 
strategies. Figure 5 sheds light over which of the two types of reasons dominate discov- 
erability for the SUSY points investigated. It shows points in the (nig, m^ghtest q) plane, 
color-coded to describe the luminosity needed for discovery in the scenario of y/s = 14 TeV 
with a combination of the MET based and the optimized Mt2 based search strategy and 
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Energy 


7 TeV 


14 TeV 


Background uncertainty cr^/b = 


0% 


10% 


50% 


0% 


10% 


50% 


C = 35 pbr 1 


■ Lv - L l z 


0.3% 


0.2% 


0.0% 


6.6% 


4.6% 


0.7% 


MET 


0.0% 


0.0% 


0.0% 


5.3% 


5.3% 


5.3% 


combined 


0.3% 


0.2% 


0.0% 


7.4% 


6.2% 


5.3% 


C = 100 ptr 1 


JW J 2 


0.7% 


0.5% 


0.1% 


16.3% 


12.5% 


3.9% 


MET 


0.0% 


0.0% 


0.0% 


16.7% 


16.7% 


16.7% 


combined 


0.7% 


0.5% 


0.1% 


19.9% 


17.6% 


16.8% 


C = 1 fbr 1 


JW J 2 


7.2% 


4.1% 


0.8% 


54.3% 


47.6% 


29.7% 


MET 


3.9% 


3.9% 


3.9% 


75.2% 


75.2% 


75.2% 


combined 


7.7% 


5.4% 


4.1% 


76.2% 


76.2% 


75.4% 


C = 10 ftr 1 


Mro 


29.9% 


26.1% 


19.8% 


85.8% 


81.6% 


65.2% 


MET 


38.4% 


38.4% 


38.4% 


95.9% 


95.9% 


95.9% 


combined 


42.5% 


42.5% 


41.6% 


96.1% 


96.0% 


96.0% 


£ = 100 fb" 1 


M T2 


(65.0%) 


(64.6%) 


(56.2%) 


96.2% 


92.7% 


85.4% 


MET 


(78.8%) 


(78.8%) 


(78.8%) 


99.5% 


99.5% 


99.5% 


combined 


(80.3%) 


(80.3%) 


(80.0%) 


99.6% 


99.6% 


99.6% 


C = 1 ab- 1 


M T2 


(82.6%) 


(80.3%) 


(72.6%) 


99.0% 


95.9% 


87.2% 


MET 


(91.6%) 


(91.6%) 


(91.6%) 


99.8% 


99.8% 


99.8% 


combined 


(92.2%) 


(92.2%) 


(92.1%) 


99.9% 


99.9% 


99.9% 



Table 3: Percentage of discoverable points (out of 1941) for the optimised Mt2 based, the 
MET based and the combined search strategies, different luminosities and y/s = 7, 14 TeV. 
The systematic uncertainty in the background is taken to be a fixed percentage, / = 0%, 
10% or 50% of the expected number of background events b, but at least one event, i.e. 
<j~ b = max |/ x b events, 1 event j. The entries in brackets show the fractions of discoverable 
points for luminosities which we do not expect to be recorded at those energies. The number 
of points discovered with our MET based strategy does not depend on the systematic 
uncertainty in the background due to none of the simulated background events passing the 
stringent cuts of this analysis. A more careful analysis [8, 9] of the expected backgrounds 
shows a strong dependence of discoverability on the size of the systematic uncertainty. 

ignoring systematic uncertainties in the background. The value Mj^ = 320 GeV opti- 
mizes the cut "in the long run", i.e. for 1 fb _1 and above. The points which are hard to 
discover predominantly populate high- mass regions of the (nig, miightest q) plane. The tran- 
sition from black (easily discoverable with ~ 1 pb" 1 ) in the bottom left corner to yellow 
(hard to discover, need > 10 fb _1 ) at the extremities does not seem to be single-valued 
and figure 6 examines this further. It shows the luminosity needed for discovery in the 
(Am, miightest sparton) plane, with Am being the difference between the mass of the lightest 
sparton and that of the LSP. Figure 6 shows that the heavier the sparticles, the harder 
it is to discover the point the sparticles belong to. However, for the same mass value for 

c light gray in a b/w printout 
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Figure 5: log 10 [luminosity (pb) _1 needed for discovery] in the (ra§ , mughtcst q) plane at 
14 TeV with a combination of the MET based and optimised Mti based strategies. The Mt2 
based strategy was optimized for an integrated luminosity of 1 fb -1 . Systematic uncertainties 
in the background have been neglected. 



the lightest sparton on the abscissa, some points are still much harder to discover than 
others. The variation of luminosity with Am suggests that another important factor for 
discoverability is the mass degeneracy of the lightest sparton and the LSP. 

A few representative spectra of hard-to-discover points are shown in figure 7, as plotted 
with PySLHA [23]. We identified 57 points which require more than 10 ftr 1 at y/a = 14 TeV 
to be discovered. Several categories can be distinguished: 

1) Points which due to high sparton masses have a small sparton production cross- 
section of less than 10 fb, while the overall SUSY cross-section is above 10 fb due to slepton 
and gaugino production. These points cannot be easily discovered at the LHC with any of 
the ATLAS MET search strategies unless dedicated studies [24, 25] are employed. These 
points might be easier to discover at a future linear collider. A representative mass spectrum 
of the 16 points identified in this category is shown in figure 7a. 

2) Points which due to the high sparticle masses have a small overall production cross- 
section of less than 10 fb. These points might prove difficult to discover at a future linear 
collider as well. A representative mass spectrum of the 15 points that fall in this category 
is shown in figure 7b. 

3) Points which are not cross-section limited, i.e. a > 10 fb. 
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Mass of lightest sparton (GeV) 

Figure 6: log 10 [luminosity (pb) _1 needed for discovery] with the combined optimal Mt2 and 
MET based strategy at t/s = 14 TeV in the (Am, m^ghtest sparton) plane. Am is the mass 
difference between the lightest sparton and the LSP. The Mti based strategy was optimized 
for an integrated luminosity of 1 fb _1 . Systematic uncertainties in the background have been 
neglected. 



The above three categories are mutually exclusive. However, the points can also be 
categorised depending on the mass degeneracies of the lightest sparton and the LSP, and 
we show such spectra in figures 7c, 7d. Any other squarks and the gluino are far more 
massive and their production is therefore suppressed. Even if the production cross-section 
for the nearly mass-degenerate squarks is large enough, above 10 fb, the mass degeneracy 
leads to a decay topology which only gives soft jets, making the signal difficult to distinguish 
from background. We identified in our sample 16 points which have a lightest squark of 
the first two generations nearly mass-degenerate with the LSP and 7 points which have a 
b or t squark nearly mass-degenerate with the LSP. 6 points fall in both categories, having 
both the lightest squark of the first two generations and the lightest squark of the third 
generation nearly mass-degenerate with the LSP. Points with mass near-degenerate particles 
could be dealt with by looking at events where the produced sparticle pair is boosted by 
the emission of hard-QCD jets from initial or final state radiation [26, 27]. 

Five hard-to-discover points are both cross-section limited and exhibit large mass de- 
generacies, which will make them particularly hard to tackle at the LHC. 
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(a) An example of a point where all squarks and 
the gluino are very heavy. Gaugino production 
dominates. (16 such points identified). 
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(c) An example of a point with nearly 
mass-degenerate squark of the first two 
generations and LSP (16 points identified). 
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(b) An example of a point where all sparticles are 
very heavy (15 such points identified). 
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(d) An example of a point with nearly 
mass-degenerate squark of the third generation 
and LSP (7 points identified). 



Figure 7: Some hard to discover points. The dotted lines show possible decay modes. Only 
branching ratios of above 5% have been drawn. The darker shades correspond to higher 
branching ratios. 



11 Discussion 

Recent papers by Conley et al. [8, 9] describe similar studies to the present one, applying 
several other ATLAS MET based search strategies on a set of ~71k model points in pMSSM. 
The similarities and differences between these studies and ours are discussed in this section. 

In the scenarios studied by Conley et al., an additional particle outside the MSSM 
is implicitly assumed to make up part or all of the dark matter relic density. This has 
the consequence of a preponderance of wino dominated lightest supersymmetric particles, 
leading to quasi-mass degenerate lightest charginos and lightest neutralinos. Reference [9] 
did not perform a global fit, but instead scanned for points which were compatible with 
each experimental 95% CL. Depending on the scan prior used, the sparticle masses were 
constrained to be less than 3 TeV or 1 TeV. Each point was given equal weight in their 
results. 
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In comparison to our work, in the MET based strategies employed by Conley et al., 
the lepton multiplicity required is varied. However, it is suggested both in their work and 
in ATLAS studies [1, 3] that the lepton channel is the most robust of the MET-based 
analyses. The backgrounds in references [8, 9] were produced by the ATLAS SUSY group 
with state-of-the-art Monte Carlo event generators and the full ATLAS detector simulation. 
The signal was generated using a fast detector simulation. 

For different scenarios with systematic uncertainty in the background of 50% or 20%, 
Conley et al. calculate how many of the generated SUSY points are discoverable using any 
of the ATLAS 11 MET analyses, requiring the same criteria for discoverability as the ones 
in the present paper. Conley et. al found an optimal cut on the value of the M c g variable 
employed in MET type analyses, while we did the same for the Mt2 variable. It is difficult 
to perform a direct analysis between our results and those by Conley et al., since their paper 
starts out from different assumptions. We note for example that setting the sparticle mass 
upper range to 1 or 3 TeV would have further increased the fraction of points discovered in 
our study. 

Reference [9] found that, at y/s = 7 TeV and after 1 fb _1 , MET cuts would discover 
between 50-90% of the models considered, depending on the systematic error assumed 
on the background. Here, we found a similar level of dependence of M^-based search 
discoverability on the systematic error assumed. References [8, 9] did not consider the Mt2 
based search strategy. 

Conley et al. also study the effect of a detector-stable charged LSP, whereas the 
LSP in the points we investigated was Xi throughout. We did find though 7 points with 
pseudo-stable xt> 1 point with pseudo-stable X2 ana ^ 13 points with pseudo-stable fi. d 
These points exhibit extreme mass degeneracies between the pseudo-stable, next-to-lightest 
supersymmetric particle and the LSP. The mass difference is a few hundred MeV and the 
decay is strongly phase-space suppressed. The pseudo-stable particles have decay lengths 
from a few centimeters to several meters. Whilst the simulation that we performed becomes 
invalid for these points since it includes no information on the detector geometry, the high 
production cross-section for these points combined with a search for long-lived charged 
particles would make such points detectable by other means [1]. 

The analysis of the discovery reach of a search strategy based on Mt2 is unique to 
our paper. Also, in our study we have taken into consideration a bigger range of sparticle 
masses, allowing for much higher hierarchies in the points studied. The limited range of 
sparticle masses restricts the type of hard-to-discover points found by references [8, 9]. Both 
our study and references [8, 9] investigated the effects of varying the systematic uncertainty 
in the background and the luminosity, but Conley et al. did not explore the discovery 
potential at yfs = 14 TeV for luminosities higher than 1 fb , a luminosity which will 
certainly be surpassed at the LHC. 



d We consider a particle to be "pseudo-stable" if it has a lifetime larger than 100 picoseconds or equivalent ly 
its decay length is longer than 3 cm. 
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12 Conclusions 



We have simulated searches for SUSY at the LHC, assuming different luminosities and 
centre of mass energies. We used a selection based on the Mt2 variable alone, an ATLAS 
MET-type search strategy and a combined approach and applied them to ~2000 pMSSM 
good-fit points. Based on these we have suggested optimal values for the cut on M?2 for 
low luminosities, which are 320 GeV for yfs = 7 TeV and 400 GeV for y/s = 14 TeV. 

Ignoring systematic uncertainties in the background and using a combined optimised 
Mt2 an d MET based search, we reach the following conclusion. If we assign equal prob- 
ability to each model studied, a 7 TeV LHC would discover SUSY with 42.5% probability 
after accumulating 10 fb _1 . Similarly, a 14 TeV LHC would discover SUSY in 96.1% of the 
cases after an integrated luminosity of 10 fb _1 . The few points that still evade detection 
are characterised by heavy gluino masses above 2.5 TeV and heavy squark masses or high 
mass degeneracies between the lightest squark and the LSP. 

The Mt2 based strategy is better suited than the MET based strategy for integrated 
luminosities up to several inverse femtobarns, in the absence of systematic uncertainties in 
the background. The Mt2 strategy is based on a single cut on Mt2- It accepts both more 
signal and more background than the MET strategy, meaning that at higher luminosities 
it is more vulnerable to fixed fractional systematic uncertainties in the background. The 
MET based strategy is more sensitive than the Mt2 based strategy to points with high 
gluino masses. Of course, a combined strategy is the most powerful. 

An optimization of the standard ATLAS searches (i.e. finding optimal values for the 
cuts on the selection variables used therein) would be desirable, though more difficult to 
make, due to the multi-dimensionality of the parameter space. A further optimization 
only makes sense in the context of good understanding of the systematic uncertainties, an 
important concern for the experiments. 

While the LHC machine and the ATLAS and CMS experiments have performed re- 
markably well, SUSY searches at the LHC are still in their infancy. Our work suggests that 
there is still much more of the parameter space to be explored and that this can be done 
efficiently using the search strategies available. 

Acknowledgments 

We are grateful to the authors of [10] for providing us with the points analysed. Sasha 
Caron has kindly provided the significance calculator for section 8. A.D. is thankful to 
Merton College, Oxford for hospitality. This work has been partially supported by STFC. 

References 

[1] The ATLAS Collaboration, G. Aad et. al, Expected Performance of the ATLAS 
Experiment - Detector, Trigger and Physics, hep-ex/0901 . 0512. 

[2] G. L. Bayatian et. al, CMS physics Technical Design Report, Volume II: Physics 

Performance. oai:cds.cern.ch:942733, J. Phys. G 34 (2006), no. CERN-LHCC-2006-021. 
CMS-TDR-008-2 995-1579. 669 p. 



- 20 - 



[3] The ATLAS Collaboration, G. Aad et. al., Search for squarks and gluinos using final states 
with jets and missing transverse momentum with the ATLAS detector in sqrt(s) = 7 TeV 
proton-proton collisions, hep-ex/1102 . 5290. 

[4] CMS Collaboration, V. Khachatryan et. al., Search for Supersymmetry in pp Collisions at 7 
TeV in Events with Jets and Missing Transverse Energy, hep-ex/1101 . 1628. 

[5] C. G. Lester and D. J. Summers, Measuring masses of semiinvisibly decaying particles pair 
produced at hadron colliders, Phys. Lett. B463 (1999) 99-103, [hep-ph/9906349]. 

[6] A. Barr, C. Lester, and P. Stephens, m(T2) : The Truth behind the glamour, J. Phys. G29 
(2003) 2343-2363, [hep-ph/0304226]. 

[7] A. J. Barr and C. Gwenlan, The race for supersymmetry: using Mt2 for discovery, Phys. 
Rev. D80 (2009) 074007, [hep-ph/0907. 2713]. 

[8] J. A. Conley, J. S. Gainer, J. L. Hewett, M. P. Le, and T. G. Rizzo, Supersymmetry Without 
Prejudice at the LHC, hep-ph/1009 . 2539. 

[9] J. A. Conley, J. S. Gainer, J. L. Hewett, M. P. Le, and T. G. Rizzo, Supersymmetry Without 
Prejudice at the 7 TeV LHC, hep-ph/1103 . 1697. 

[10] S. S. AbdusSalam, B. C. Allanach, F. Quevedo, F. Feroz, and M. Hobson, Fitting the 
Phenomenological MSSM, Phys. Rev. D81 (2010) 095012, [hep-ph/0904 . 2548]. 

[11] F. Feroz and M. P. Hobson, Multimodal nested sampling: an efficient and robust alternative 
to MCMC methods for astronomical data analysis, astro-ph/0704.3704. 

[12] F. Feroz, M. P. Hobson, and M. Bridges, MultiNest: an efficient and robust Bayesian 
inference tool for cosmology and particle physics, astro-ph/0809 . 3437. 

[13] B. C. Allanach et. al, The Snowmass points and slopes: Benchmarks for SUSY searches, 
Eur. Phys. J. C25 (2002) 113-123, [hep-ph/0202233]. 

[14] The ATLAS Collaboration, G. Aad et. al., Prospects for Supersymmetry discovery based on 
inclusive searches at a 7 TeV centre- of-mass energy with the ATLAS detector, Tech. Rep. 
ATL-PHYS-PUB-2010-010, CERN, Geneva, Jul, 2010. 

[15] H.-C. Cheng and Z. Han, Minimal Kinematic Constraints and MT2, JEEP 12 (2008) 063, 
[hep-ph/0810.5178]. 

[16] The ATLAS Collaboration, G. Aad et. al., Early supersymmetry searches in channels with 
jets and missing transverse momentum with the ATLAS detector, Tech. Rep. 
ATLAS-CONF-2010-065, CERN, Geneva, Jul, 2010. 

[17] M. Bahr et. al., Herwig++ Physics and Manual, Eur. Phys. J. C58 (2008) 639-707, 
[arXiv: 0803. 0883]. 

[18] S. Gieseke et. al, Herwig++ 2.5 Release Note, hep-ph/1102 . 1672. 

[19] P. Z. Skands et. al., SUSY Les Houches Accord: Interfacing SUSY Spectrum Calculators, 
Decay Packages, and Event Generators, JEEP 07 (2004) 036, [hep-ph/0311123]. 

[20] M. Cacciari and G. P. Salam, Dispelling the N 3 myth for the k t jet-finder, Phys. Lett. B641 
(2006) 57-61, [hep-ph/0512210]. 

[21] M. Cacciari, G. P. Salam, and G. Soyez, The anti-kt jet clustering algorithm, JEEP 04 
(2008) 063, [hep-ph/0802.1189]. 

[22] R. D. Cousins, J. T. Linnemann, and J. Tucker, Evaluation of three methods for calculating 



- 21 - 



statistical significance when incorporating a systematic uncertainty into a test of the 
background- only hypothesis for a Poisson process, Nucl. Instrum. Meth. A595 (2008) 
480-501. 



[23] A.Buckley, http://www.insectnation.org/projects/pyslha. 

[24] E. Lytken, Searches for Direct Slepton Production with ATLAS, Czech. J '.Phys. 54 (2004) 
A169-A173. 

[25] C.-Y. Chen and A. Freitas, General analysis of signals with two leptons and missing energy 
at the Large Hadron Collider, JHEP 1102 (2011) 002, [hep-ph/1011 . 5276]. 

[26] D. S. M. Alves, E. Izaguirre, and J. G. Wacker, Where the Sidewalk Ends: Jets and Missing 
Energy Search Strategies for the 7 TeV LEG, hep-ph/1102 . 5338. 

[27] J. Alwall. M.-P. Le, M. Lisanti. and J. G. Wacker, Searching for Directly Decaying Gluinos at 
the Tevatron, Phys. Lett. B666 (2008) 34-37, [hep-ph/0803 . 0019]. 



- 22 - 



